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Abstract 

Husimi distribution function for the one-dimensional Ising model is 
obtained. One-point and joint distribution functions are calculated 
and their thermal behaviour are discussed. 
PACS numbers: 02.70.-c, 03. 65. Ca 

1 Introduction 

Husimi distribution function P , is a positive definite distribution over the pa- 
rameter space of a definite set of coherent states 2]. Positive definitness of the 
Husimi distribution function, makes it as a good candidate for investigating 
quantum features using a classical hke distribution over the phase space. In 
this case, to each coherent state \z) of the system there corresponds a unique 
distribution function = ^{zle'^^^lz), where Z, is the partition function 
of the system. H and /3 are Hamiltonian and inverse temperature parameter 
respectively. In this paper, using the single-fermion or spin-| coherent states, 
the coherent states of a one-dimensional spin chain are obtained. The Husimi 
distribution function is then obtained for the one-dimensional Ising model as 
the prototype of some important models in statistical mechanics. One-point 
and joint distribution functions or correlation functions are calculated and 
their thermal behaviour is discussed. 
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2 One-dimensional Ising model 

Hamiltonian of the one dimensional Ising model in a homogeneous magnètic 
field B along the ^-axis is 

k=l k=l 

where periòdic boundary condition is assumed, i.e, + 1 = 1. The Hilbert 
space of this model is a 2^-dimensional vector space spaned by the foUowing 
tensor products as basis vectors 

N 

\i) = |zi) ^ \Í2) ® ... ^ \Ín) = n ^h), (2) 

k=l 

where \ik) G {|+), |— )}, for k ~ 1...N. The basis vectors are Hamiltonian 
eigenvectors with respective eigenvalues -E^, 

N N 

H\i) = {-J^ikik+i + B^ik)\i), 

k=l k=l 
N N 

= -J^ik^k+i^ B^ik, (3) 

k=l k=l 

3 Ising model coherent states 

Coherent states of the one dimensional Ising model can bc dcfined from the 
tensor product of the single-fermion coherent states [2]. A single-fermion 
coherent state for the spin half {s — ^) representation can be written like 
this 

9 9 
1^) =sin(-)exp(-z(/p)|+)+cos(-)|-), (4) 

where z = {^)exp{—i(f), 0<9<7r, 0<(^<2tt and |+), |— ) are Sz 
eigenvectors. The completeness of the fermion coherent states is 

^ ^dü\z){z\ = l, (5) 



2tt 

where dfl — sin 9 d9d(p. Tensor product of these fermion coherent states 
makes Ising model coherent states 

Na a 

|i) = n®(sin(^)expHv..)|+) + cos(^)|-)), (6) 
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the scalar product between a coherent state \z> and an eigenvector \i >, of 
the Hamiltonian H, is 

{z\i) = exp(-i(^fe)5í„i + cos(y)(Jí, _i), (7) 

so 

N 



= n(sm^(y)^^.,i + cos^íy)^..,-!)- (8) 



Partition function of the Hamiltonian (1) is defined as 

Z = ^(í I eM-m I ^) = íKÍ''^), (9) 



where 



^ _ íeMí^iJ - B)) exp(-/3J) \ 

exp(-/3J) exp(/3(J+S))i ' ^^""^ 



is the transition matrix. 



4 Husimi distribution function 

Husimi distribution function is a normahzed positive definite distribution in 
the parameter space of coherent states and is defined as foUows 

H{z) = ^{z\exp{-pH)\z>, 

= ^E(^exp(-/?^)|ï)(ïlf>, 

= l^e^·^S^=i'^'^+i-^^Sr=i^fc|(i^ï)|2, (11) 
^{?} 

substituting (8) in (11) and doing some routine calculations, we get the 
foUowing relation for Husimi distribution 

1 ^ . . . 

m=l m<n 
+ {-l)^{Si^zS2,z···SN,z)UlU2···Un}, (12) 
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where for simplicity = cosO^ is assumed and a general n-point corelation 
function {Sm,zSn,z • ■ • Sr,z) is defined as 

{Sm,zSn,z---Sr,z) = ^^t^t^ ■ ■ ■ trC^'^L·^^'^^^-^^EL·^^ , (13) 

i 

where Sk^zlik) = iklik), from (12) it is clear that the Husimi distribution 
function is independent of parameters (pk and contains the information about 
all the n-point correlation functions, i.e., it is the generator of the n-point 
functions of the Ising model 

{Sri,z^r2,z ■ ■ ■ Srm,z) — (—3)"* ' ' ' l·^{^)'^ri'^r2 ' ' ' "^rm H ^"^kj 

k=l 

= (-3r{Ur,Ur,---UrJ. (14) 

Now the Husimi distribution function rclatcd to the cite k, can be calculated 
by integrating over parameters of the phase spase except those belong to the 
kth cite, so 

/^K) = 2(i-«.^[:j^), (15) 



where as usual is 



= J ■ (16) 



The joint Husimi distribution function related to cites m and n, is equal to 

, , 1, , Ar[I'^a,] tr[f'^-"+"'a-f"-"'a-]^ , ^ 

^{Um,Un) = -{l-{Um + Un)-^^^l^;^ + UmUn ^^^j^;^^^^ ), (17) 

for calculating ír[T(T^] and tr\T^~^'^'^azT'^~'^o'z\ we can diagonalize the ma- 
trix T, and write the matrix cr^ in the basis of eigenvectors of T, the eigen- 
values of T are 



A± = exp(/?J) cosh(/3S) ± ·\/exp(2/3J) cosh.\(3B) - 2 sinh(/?J), (18) 

and the transformation matrix between the Standard basis {|+),|— )} and 
the eigenvectors of T is 

t/=HÍ -cos(^)\ 
\cos(a;j sm(a;) / 



where |A+) = í ■^^'^^.'^-^ j g^^^j |^_\ _ í are normalized eigenvectors 

of T belonging to eigenvalues A+ and A_ respectively and 

tanM = /"P^-^^^ (20) 

exp(/3 J) sinh(/3B) + -^exp(2/?J) smh^(/35) + exp(-2/3J) 

Now in the basis |A+), |A_), we have 

trlfa] - íríK W- cos(2a;) -sin(2c.)V 
tr[la4 - tr[y^ A^J l^-sin(2u;) cos(2a;)j^' 

J-X^cos{2üj) -A^sin(2a;)Y 
~ *'"^l^-A^sin(2a;) A^cos(2a;) j^' 

= (A^- A^)cos(2a;), (21) 



and similarly 



ír[f^-^'+V,f^-V,] = A^cos'(2cu) + AÍ'--''+*Ar'sin2(2L^) 

+ A^-^+*A^+-*sin2(2a;) + A^cos2(2a;), (22) 

substituting (21) and (22) in (15) and (17) respectively, we obtain 



1 - 
fx{uk) = -(1 - cos(2a;)[ ^^^^+ ]^fc), (23) 



1 A^ - A^ 

n{ui,Uj) = -{1- {ui + Uj)cos{2uj)[^^j-^-^]uk + UiUjCOS^{2uj) 

, \N-j+^^.j-^ 

+ UiUjSin'{2uj)^ AF+Af ^ ^ 

in thermodynamic limit, i.e., A*" — > oo, we have 

/u(Mfc) = ^(1 + COs(2w)Mfc), 

li{ui, Uj) — -(1 + (tíj + Uj) cos(2ü;) + UiUj cos^(2a;) 

+ «,«,(^)^-^sin2(2a;)), (25) 



such that 



cos(2a;) 



2 exp(/5 J) smh{pB) 



^4exp(2/9J) sinh2(/3S) + 4exp(-2/3J) ' 




exp(-4/3J) ' 

sinh^(/3_B) 



tan(2a;) 



exp{-2(3J) 
smh(/3B) 



(26) 



where Sgn{B), is the sign function. The behaviour of one-point distribution 
function ii{uk), in high and low temperatures, is 



in zero tempcraturc, //(m^) attains it's maximum for — 0, which from (4) 
corresponds to the spin down state, |— ), as expected. 

Similarly we can find the foUowing behaviour for joint distribution //(wj, Uj), 



in this case the disjoint distribution attains it's maximum value for 9i — 6j — 

0, which means that most probably the two spins are in the spin down state, 

1. e., |— ), in derivation of these behaviours it is assumed that Sign[B) — 1, 
by reversing the magnètic field {Sign{B) = — 1), the spins flip to the spin up 
state. 
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